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Abstract 

A single vortex in the charged Bose gas (CBG) has a charged core and its profile 
different from the vortex in neutral and BCS superfluids. Lower and upper critical 
fields of CBG are discussed. The unusual resistive upper critical field, H C 2(T), of 
many cuprates and a few other unconventional superconductors is described as the 
Bose-Einstein condensation field of preformed bosons-bipolarons. Its nonlinear tem- 
perature dependence follows from the scaling arguments. Exceeding the Pauli para- 
magnetic limit is explained. Controversy in the determination of H C 2(T) of cuprates 
from kinetic and thermodynamic measurements is addressed in the framework of 
the bipolaron theory. 
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Introduction 

The seminal work by Bardeen, 
Cooper and Schrieffer [1] taken fur- 
ther by Eliashberg [2] to the interme- 
diate coupling solved one of the ma- 
jor problems in Condensed Matter 
Physics. High-temperature super- 
conductors present a challenge to the 
conventional theory. While the BCS 
theory provides a qualitatively cor- 
rect description of some novel super- 
conductors like magnesium diborade 
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and doped fullerenes (if the phonon 
dressing of carriers, i.e. polaron for- 
mation is properly taken into ac- 
count), cuprates remain a problem. 
Here strong antiferromagnetic and 
charge fluctuations and the Frohlich 
and Jahn- Teller electron-phonon in- 
teractions have been identified as an 
essential piece of physics. In partic- 
ular, experimental [3,4,5,6,7,8,9,10] 
evidence for an exceptionally strong 
electron-phonon interaction in all 
high temperature superconductors is 
now overwhelming. Our view, which 
we discussed in detail elsewhere [11] 
is that the extension of the BCS the- 
ory towards the strong interaction 
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between electrons and ion vibrations 
describes the phenomenon naturally. 
The high temperature superconduc- 
tivity exists in the crossover region 
of the electron-phonon interaction 
strength from the BCS-like to bipo- 
laronic superconductivity as was 
predicted before [12], and explored 
in greater detail after the discovery 
[13,14,15,16]. The low energy physics 
in this regime is that of a charged 
Bose gas of small bipolarons, which 
are real-space bosons dressed by 
phonons. They are itinerant quasi- 
particles existing in the Bloch states 
at temperatures below the character- 
istic phonon frequency. Here I review 
the bipolaron theory of the vortex 
state. 



1 Charged vortex 



CBG is an extreme type II supercon- 
ductor, as shown below. We can anal- 
yse a single vortex in CBG and calcu- 
late the critical fields by solving a sta- 
tionary equation for the macroscopic 
condensate wave function if) s (r) [17], 
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Subtracting rib m the integral of 
Eq.(l) explicitly takes into account 
the Coulomb interaction with the 
homogeneous charge background 
of the same density as the density 
of charged bosons rif,. Here 2e and 
m** are the charge per boson and 
the effective mass, respectively, and 



h — c — ks — 1). 

The integra-differential equation (1) 
is quite different from the Ginsburg- 
Landau [18] and Gross-Pitaevskii 
[19] equations, describing the vortex 
in the BCS and netral superfluids, 
respectively. While CBG shares the 
quantum coherence with the BCS 
superconductors and neutral super- 
fluids owing to the Bose-Einstein 
condensate (BEC), the long-range 
(nonlocal) interaction leads to some 
peculiarities. In particular, the vor- 
tex is charged in CBG, and the co- 
herence length is just the same as 
the screening radius. 



Indeed, introducing dimensionless 
quantities / = \ip s \/nj 2 , p = r/A(0), 
and h =2e£(0)A(0)V x A for the or- 
der parameter, length and magnetic 
field, respectively, Eq.(l) and the 
Maxwell equations take the following 
form: 
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A new feature compared with the 
GL equations for a single vortex [20] 
is the electric field potential deter- 
mined as 



(b=-^— I dv'Vir-r') 
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with a new fundamental unit <p c = 



2 



em**£(0) 2 . The potential is calcu- 
lated using the Poisson equation (3). 
At T = the coherence length is the 
same as the screening radius, 



f (0) = {2 1 ' 2 m**u ps )- 1 '\ 



(6) 



and the London penetration depth is 

(7) 
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Here u ps = [16ne 2 n b / (e m**)] 1/2 is 
the plasma frequency. There are now 
six boundary conditions in a single- 
vortex problem. Four of them are the 
same as in the BCS superconductor 
[20], h = dh/p = 0, / = 1 for p = oo 
and the flux quantization condition, 
dh/dp = —pf 2 /np for p — 0, where p 
is an integer. The remaining two con- 
ditions are derived from the global 
charge neutrality, = for p = oo 
and 



0(0)= p\n{p){l-f)dp 



stant the Coulomb repulsion remains 
weak even for heavy bipolarons, 

4m**e 2 . . 

If k > 1, Eq.(4) is reduced to the 
London equation with the familiar 
solution h = pK (p)/n, where K (p) 
is the Hankel function of imaginary 
argument of zero order. For the re- 
gion p < p, where the order param- 
eter and the electric field differ from 
unity and zero, respectively, we can 
use the flux quantization condition 
to "integrate out" the magnetic field 
in Eq.(2). That leaves us with two 
parameter-free equations written for 
r = Kp as 



1 d df p 2 f 

~i~ r i r _ vf 

r dr dr r z 
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and 

1 d d(f) 
r dr dr 
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for the electric field at the origin, p = 

0. We notice that the chemical po- 
tential /i is zero at any point in the 
thermal equilibrium. 

CBG is an extreme type II supercon- 
ductor with a very large Ginsburg- 
Landau parameter, k = A(0)/£(0) 3> 

1. Indeed, with the material parame- 
ters typical for oxides, such as m** = 
10m e , rib = 10 21 cm~ 3 and the static 
dielectric constant eo = 10 3 we ob- 
tain f (0) ~ 0.48nm, A(0) ~ 265nm, 
and the Ginsburg-Landau ratio k ~ 
552. Owing to a large dielectric con- 



They are satisfied by regular so- 
lutions of the form / = c p r p and 
= 0(0) + (r 2 /4), when r -> 0. The 
constants c p and 0(0) are determined 
by complete numerical integration 
of Eqs.(10) and (11). The numerical 
results for p = 1 are c\ ~ 1.5188 and 
0(0) ~ -1.0515. 

In the region p « r < pn the so- 
lutions are f — 1 + (Ap 2 /r 4 ) and 
= —p 2 /r 2 . In this region / dif- 
fers qualitatively from the BCS or- 
der parameter, fscs — 1 — (p 2 /r 2 ) 
[20]. The difference is due to a local 
charge redistribution caused by the 
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magnetic field in CBG. Quite differ- 
ent from the BCS superconductor, 
where the total density of electrons 
remains constant across the sample, 
CBG allows for flux penetration by 
redistributing the density of bosons 
within the coherence volume. This 
leads to an increase of the order pa- 
rameter compared with the homoge- 
neous case (/ = 1) in the region close 
to the vortex core. Inside the core 
the order parameter is suppressed, as 
in the BCS superconductor. The re- 
sulting electric field, (together with 
the magnetic field) acts as an addi- 
tional centrifugal force increasing the 
steepness (c p ) of the order parameter 
compared with the BCS superfluid, 
where C\ ~ 1.1664. 

The breakdown of the local charge 
neutrality is due to the absence of 
any equilibrium normal state solu- 
tion in CBG below H c2 {T) line. Both 
superconducting (A k 7^ 0) and nor- 
mal (Ak = 0) solutions are allowed 
at any temperature in the BCS su- 
perconductors. Then the system de- 
cides which of two phases (or their 
mixture) is energetically favorable, 
but the local charge neutrality is re- 
spected. In contrast, there is no equi- 
librium normal state solution (with 
i/j s = 0) in CBG below # c2 (T)-line 
because it does not respect the den- 
sity sum rule. Hence, there are no 
different phases to mix, and the only 
way to acquire a flux in the thermal 
equilibrium is to redistribute the lo- 
cal density of bosons at the expense 
of their Coulomb energy. This energy 
determines the vortex free energy 
F — E v — Eq, which is the difference 
of the energy of CBG with, E v , and 
without, E , magnetic flux, 



F = J dr— |[V + 2*eA(r)]^(r)| 2 
+ e0 c 0[|^(r)| 2 -n b ]+ (V g 7r A)2 . 

Using Eqs.(2), (3) and (4) it can be 
written in the dimensionless form as 

00 

F = 2irJ[h 2 - -0(1 + f)]pdp. (12) 


In the large k limit the main con- 
tribution comes from the region 
p/n < p < p, where / ~ 1 and 
(j) ~ —p 2 /( y K 2 p 2 ). The energy is thus 
the same as that in the BCS super- 
conductor, F ~ 2Tcp 2 ln(/c) / k 2 . It is 
seen that the most stable solution is 
the formation of the vortex with one 
flux quantum, p — 1, and the lower 
critical field is the same as in the BCS 
superconductor, h c \ rs ln/c/(2/c) [20]. 
However, different from the BCS su- 
perconductor, where the Ginsburg- 
Landau phenomenology is micro- 
scopically justified in the temper- 
ature region close to T c , the CBG 
vortex structure is derived here in 
the low temperature region. Actu- 
ally the zero temperature solution is 
applied in a wide temperature region 
well below the Bose-Einstein con- 
densation temperature, where the 
depletion of the condensate remains 
small. The actual size of the charged 
core is about 4£. 



2 Upper critical field in the 
strong-coupling regime 

If we "switch off" the Coulomb re- 
pulsion between bosons, an ideal 
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CBG cannot be bose-condensed at fi- 
nite temperatures in a homogeneous 
magnetic field because of a one- 
dimensional particle motion at the 
lowest Landau level [21]. However, an 
interacting charged Bose-gas is con- 
densed in a field lower than a certain 
critical value H c2 (T) [22]. Collisions 
between bosons and/or with impu- 
rities and phonons make the motion 
three-dimensional, and eliminate the 
one-dimensional singularity of the 
density of states, which prevents 
BEC of the ideal gas in the filed. 
As we show below the upper criti- 
cal field of CBG differs significantly 
from H c2 (T) of BCS superconduc- 
tors. It has an unusual positive cur- 
vature near T c , H c2 {T) ~ (T c - Tf' 2 
and diverges at T — > 0, if there is no 
localisation due to a random poten- 
tial. The localization can drastically 
change the low-temperature behav- 
ior of H c2 (T), so that at high density 
of impurities a re-entry effect to the 
normal state might occur. 

In line with the conventional defini- 
tion, H c2 (T) is a field, where a first 
nonzero solution of the linearized sta- 
tionary equation for the macroscopic 
condensate wave function occurs, 



2m 
V scat (r)ip s (r) 



[V - 2ieA(r)] 2 + /i 



(13) 



in the homogeneous magnetic field is 

fa ~ 



e n = uj(n + 1/2) + 



2m** ' 



(14) 



where u = 2eH c2 /m** and n = 
0, 1,2, ...oo. BEC occurs when the 
chemical potential "touches" the 
lowest band edge from below, i.e. 
fx — cu/2. Hence, quite different from 
the GL equation, the Schrodinger 
equation (13) does not allow for a di- 
rect determination of H c2 , In fact, it 
determines the value of the chemical 
potential. Then using this value the 
upper critical field is found from the 
density sum rule, 

oo 

J f(e)N(e,H c2 )de = n b , (15) 

E c 



where N(e, H c2 ) is the density of 
states (DOS) of the Hamiltonian, 
Eq.(13),/( e ) = [exp(e-/i)/T-l]- 1 
is the Bose-Einstein distribution 
function, and E c is the lowest band 
edge. For ideal bosons we have 
/i — E c — uj/2 and 



N(e,H c2 ) 



V2(m**] 



3/2 



CO 



4tt 2 



oo 



n=o \Je — oj{n + 1/2) 



Substituting this DOS into Eq.(15) 
yields 



Here we introduce the "scattering" 
potential V scat (r) caused, for ex- 
ample, by particle-particle and/or 
particle-impurity collisions. Let us 
first discuss noninteracting bosons, 
Vscat( r ) = 0. Their energy spectrum 



\/2{m**f/ 2 uj 7 dx 



4tt 2 
= n b - n(T), 

where 



x 1 / 2 exp(x/T) - 1 
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n{T) = 



x 



y/2(m**) 3 / 2 
4tt 2 

OO 1 

n=l V- 



oo 

/ 



a; — urn 



dx 

exp(x/T) - 1 
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is the number of bosons occupying 
the levels from n = 1 to n = oo. 
This number is practically the same 
as in zero field, n(T) = nb(T/T c ) 3 ^ 2 , if 
uj <C T c . On the contrary, the number 
of bosons on the lowest level, n — 0, 
is given by a divergent integral on the 
left-hand side of Eq.(18). Hence the 
only solution to Eq.(16) is H c2 (T) = 
0. 

The scattering of bosons effectively 
removes the one-dimensional singu- 
larity in iV (e, H c2 ) oc uo(e - u/2)- 1 / 2 
leading to a finite DOS near the bot- 
tom of the lowest level, 



N (e,H c2 ) oc 



H. 



c2 



r (H c2 ) 



(18) 



Using the Fermi-Dirac golden rule 
the collision broadening of the lowest 
level T (H c2 ) is proportional to the 
same DOS 



r (# c2 )c<iVo(e,# c2 ), 



(19) 



HJT), 
arbitrary 
units 




Fig. 1. Upper critical field of CBG com- 
pared with H c2 (T) of BCS supercon- 
ductors. 

tegral at To from below. Finally we 
arrive at 

H c2 {T) = H {t- 1 -t l ' 2 fl\ (21) 



where t = T/T c , and H is a tempera- 
ture independent constant. The scal- 
ing constant H depends on the scat- 
tering mechanism. If we write H = 
$o/(27r£o), then the characteristic 
length is 




1/4 



(22) 



so that r scales with the field as 
Fo(H c2 ) oc H 2 c2 '. Then the number of 
bosons at the lowest level is estimated 



as 



n = 



A^ 2 



r 



x 1 / 2 exp(x/T) — 1 



(20) 



as long as T ^> r . Here we apply the 
one-dimensional DOS, but cut the in- 



where I is the zero-field mean-free 
path of low energy bosons. 

The upper critical field has a nonlin- 
ear behaviour, 

H c2 (T) oc (T c -T) 3 / 2 , 

in the vicinity of T c , and diverges at 
low-temperatures as 

H c2 {T) oc T~ 3/2 . 
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These simple scaling arguments are 
fully confirmed by DOS calculations 
with impurity [22] and boson-boson 
[23] scattering. The "coherence" 
length £ of CBG, Eq.(22), depends 
on the mean free path / and the 

— 1/3 

inter-particle distance n b . It has 
nothing to do with the size of the 
bipolaron, and could be as large as 
the coherence length of the weak- 
coupling BCS superconductors. 

Thus H c2 (T) of strongly-coupled su- 
perconductors has a "3/2" curvature 
near T c different from the linear BCS 
H c2 (T). The curvature is a univer- 
sal feature of CBG, which does not 
depend on a particular scattering 
mechanism and on approximations 
made. Another interesting feature 
of strongly-coupled superconductors 
is a breakdown of the Pauli param- 
agnetic limit given by H p ~ 1.84T C 
in the weak-coupling theory. H c2 {T) 
of bipolarons exceeds this limit be- 
cause the singlet bipolaron binding 
energy A is much larger than their 
T c . Bosons are condensed at T = 
no matter what their energy spec- 
trum is. Hence, in the charged Bose- 
gas model, H c2 (0) = oo, Fig.l. For 
composed bosons, like bipolarons, 
the pair-breaking limit is given by 
/i B H c2 (0) « A, so that H c2 (0) > H p . 



3 Universal upper critical field 
of unconventional supercon- 
ductors 



In cuprates [24,25,26,27,28,29,30], 
spin-ladders [31] and organic super- 
conductors [32] high magnetic field 
studies revealed a non-BCS upward 



curvature of resistive H c2 (T). When 
measurements were performed on 
low-T c unconventional superconduc- 
tors [25,26,28,31,32], the Pauli limit 
was exceeded by several times. A 
non-linear temperature dependence 
in the vicinity of T c was unam- 
bigously observed in a few samples 
[27,28,29,30]. Importantly, a thermo- 
dynamically determined H c2 turned 
out much higher than the resistive 
H c2 [33] due to contrasting magnetic 
field dependencies of the specific heat 
anomaly and of resistive transition. 

We believe that many unconven- 
tional superconductors are in the 
'bosonic' limit of preformed real- 
space bipolarons, so their resistive 
H c2 is actually a critical field of 
the Bose-Einstein condensation of 
charged bosons [22]. Calculations 
carried out for the heat capacity of 
CBG (see below) lead to the con- 
clusion that the resistive H c2 and 
the thermodynamically determined 
H c2 are very different in bosonic su- 
perconductors. While the magnetic 
field destroys the condensate of ideal 
bosons, it hardly shifts the specific 
heat anomaly as observed. 

A comprehensive scaling of resistive 
H c2 measurements in unconventional 
superconductors is shown in Fig. 2 
[30] in the framework of the mi- 
croscopic model of charged bosons 
scattered off impurities (section 2). 
Generalised Eq.(21) accounting for a 
temperature dependence of the num- 
ber of delocalised bosons, nj,(T), can 
be written as 



H c2 (T) = H 



n b {T) 
tn b (T c ) 



_ f l/2 



3/2 



.(23) 
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In the vicinity of T c one obtains 
the parameter- free H c2 (T) oc (1 — 
tf/ 2 using Eq.(23), but the low- 
temperature behaviour depends on a 
particular scattering mechanism, and 
a detailed structure of the density of 
localised states. As suggested by the 
normal state Hall measurements in 
cuprates n b (T) can be parameterised 
as n b (T) = n b (0) + constant x T 
[34], so that H c2 (T) is described by a 
single-parameter expression as 



H c2 (T) = H 



6(1-0 



+ 1 - 1 1 '* 



3/2 



The parameter b is proportional to 
the number of delocalised bosons at 
zero temperature. We expect that 
this expression is applied in the whole 
temperature region except ultra- 
low temperatures, where the Fermi 
Golden-rule in the scaling fails. Ex- 
ceeding the Pauli pair-breaking limit 
readely follows from the fact, that 
the singlet-pair binding energy is 
related to the normal-state pseudo- 
gap temperature T*, rather than to 
T c . T* is higher than T c in bosonic 
superconductors, and cuprates. 



X 



.(24) 




T/T„ 



Fig. 2. Resisitive upper critical field 
(determined at 50% of the transi- 
tion) of cuprates, spin-ladders and 
organic superconductors scaled ac- 
cording to Eq.(24). The parameter b 
is 1 (solid line), 0.02 (dashed-dotted 
line), 0.0012 (dotted line), and 
(dashed line). The inset shows a uni- 
versal scaling of the same data near 
T c on the logarithmic scale. Sym- 
bols correspond to Tl - 2201 (•), 
LSCO{A), Bi-2201(x), Bi-2212(*), 
YBCO{o), La 2 - x Ce x Cu0 4 ~y{a), 
Sr 2 Cai 2 Cu 2 404i(+), and 
(TMTSF) 2 PF 6 (V) 



The universal scaling of H c2 near T c 
is confirmed by resistive measure- 
ments of the upper critical field of 
many cuprates, spin-ladders, and 
organic superconductors, as shown 
in Fig. 2. All measurements reveal 
a universal (1 — t) 3 / 2 behaviour in 
a wide temperature region (inset), 
when they are fitted by Eq.(24). 
The low-temperature behaviour of 
H C 2(T) / H is not universal, but well 
described using the same equation 
with the single fitting parameter, b. 



The parameter is close to 1 in high 
quality cuprates with a very narrow 
resistive transition [29]. It naturally 
becomes rather small in overdoped 
cuprates where randomness is more 
essential, so almost all bosons are lo- 
calised (at least in one dimension) at 
zero temperature. 
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4 Specific heat anomaly in 
CBG 



Bose liquids (or more precisely He 4 ) 
show the characteristic A -point sin- 
gularity of their specific heat, but 
superfluid Fermi liquids like BCS su- 
perconductors exhibit a sharp second 
order phase transition accompanied 
by a finite jump in the specific heat. 
It was established beyond doubt 
[35,36,37,38,39] that the anomaly in 
high T c cuprates differs qualitatively 
from the BSC prediction. As was 
stressed by Salamon et al.[40] the 
heat capacity is logarithmic near the 
transition, and consequently, cannot 
be adequately treated by the mean- 
field BCS theory even including the 
gaussian fluctuations. In particular, 
estimates using the gaussian fluctua- 
tions yield an unusually small coher- 
ence volume [36], and Gi number of 
the order of one. 

The magnetic field dependence of 
the anomaly [41] is also unusual, but 
it can be described by the bipolaron 
model [42,30]. Calculations of the 
specific heat of charged bosons in a 
magnetic field require an analytical 
DOS, N(e, B) of a particle, scattered 
by other particles and/or by a ran- 
dom potential of impurities. We can 
use DOS in the magnetic filed with 
an impurity scattering, which allows 
for an analytical result [30]. The spe- 
cific heat coefficient 



C(T, B) d 



T 



de 



N(e,B)e 



TdTJ exp[(e- n)/T] - 




1.0 1.2 1.4 



calculated with this DOS and with ix 
determined from rib = I deN(e, B)f(e) 



Fig. 3. Temperature dependence of the 
specific heat devided by temperature 
(arb. units) of the charged Bose-gas 
scattered off impurities for several fields 
{uj h = 2eB/m**). Fig. 3b shows two 
anomalies, the lowest one traces re- 
sistive transition, while the highest 
anomaly is the normal state feature. 

is shown in Fig. 3. 

The broad maximum at T ~ T c 
is practically the same as in the 
ideal Bose gas without scattering 
[42]. It barely shifts in the mag- 
netic field. However, there is the sec- 
ond anomaly at lower temperatures, 
which is absent in the ideal gas. It 
shifts with the magnetic field, trac- 
ing precisely the resistive transition, 
as clearly seen from the difference 
between the specific heat in the field 
and zero-field curve, Fig. 3b. The 
specific heat, Fig. 3, is in striking re- 
semblance with the Geneva group's 
experiments on DyBa 2 Cu 3 7 and 
on Y Ba 2 Cu^Oi [41], where both 
anomalies were observed. Within the 
bipolaron model, when the magnetic 
field is applied, it hardly changes 
the temperature dependence of the 
j chemical potential near the zero field 
T c because the energy spectrum of 
thermally excited bosons is practi- 
cally unchanged. That is because 
their characteristic energy (of the 
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order of T c ) remains huge compared 
with the magnetic energy of the or- 
der of 2eB/m**. In contrast, the en- 
ergy spectrum of low energy bosons 
is strongly perturbed even by a weak 
magnetic field. As a result the chemi- 
cal potential 'touches' the band edge 
at lower temperatures, while having 
almost the same 'kink'-like temper- 
ature dependence around T c as in 
zero field. While the lower anomaly 
corresponds to the true long-range 
order due to the Bose-Einstein con- 
densation, the higher one is just a 
'memory' about the zero-field transi- 
tion. This microscopic consideration 
shows that a genuine phase transi- 
tion into a superconducting state is 
related to resistive transition and 
to the lower specific heat anomaly, 
while the broad higher anomaly 
is the normal state feature of the 
bosonic system in the external mag- 
netic field. Different from the BCS 
superconductor these two anoma- 
lies are well separated in the bosonic 
superconductor at any field but zero. 

In conclusion, the bipolaron theory 
of the critical fields and vortex struc- 
tures in strong-coupling supercon- 
ductors has been reviewed. A single 
vortex in this regime has a charged 
core and its profile different from the 
vortex in neutral and BCS superflu- 
ids. The upper critical field is also 
qualitatively different from the weak 
and intermediate-coupling (T) . 
We have interpreted unusual resis- 
tive upper critical fields of many un- 
conventional superconductors as the 
Bose-Einstein condensation field of 
preformed bosons-bipolarons. Their 
nonlinear temperature dependences 
follow from the scaling arguments. 



Exceeding the Pauli paramagnetic 
limit has been explained, and the 
controversy in the determination of 
H C 2(T) of cuprates from kinetic and 
thermodynamic measurements has 
been addressed in the framework of 
the bipolaron theory. 
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